Riemann surfaces, separation of variables and classical 
and quantum integrability. 



0. Babelon 1 and M. Talon Q 
Laboratoire de Physique Theorique et Hautes Energies. 



Abstract 

We show that Riemann surfaces, and separated variables immediately provide classical Pois- 
son commuting Hamiltonians. We show that Baxter's equations for separated variables immedi- 
ately provide quantum commuting Hamiltonians. The construction is simple, general, and does 
not rely on the Yang-Baxter equation. 
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1 Introduction. 



We know since Liouville that integrability means commuting Hamiltonians. It is the primary role 
of Lax matrices and the Yang-Baxter equation to provide non trivial such Hamiltonians. In the 
classical theory, additional benefits are the spectral curve T and the ability to separate variables 
by considering g = genus(r) points on it 

In the quantum theory, the analog construction is Sklyanin's method of separation of variables and 
Baxter's equations |2|, ||. Despite the beauty of this result, the route from a Yang-Baxter defined 
quantum integrable model to the separated variables is usually long and difficult, especially in the 
non hyper elliptic case. 

Here, we show that we can reverse the strategy. We start from separated variables and consider 
Baxter's equations as equations for the Hamiltonians. We then prove that these Hamiltonians 
commute under very general hypothesis. 

By its generality, its simplicity and its close analogy to the classical case, this result could provide 
a good starting point to build a theory of quantum integrable systems. 



2 The main theorem. 



Consider a curve in C 2 

a 

r(\, f i) = R (\,n) + Y,Rj(\v)H j = o (1) 

3=1 

where the Hi are the only dynamical moduli, so that Ro(X, fi) and Ri(\, h) do not contain any 
dynamical variables. If things are set up so that T is of genus g and there are exactly g Hamiltonian 
Hj (see below for realizations of this setup), then the curve is completely determined by requiring 
that it passes through g points (Aj,/Xj), i = 1, Indeed, the moduli Hj are determined by 

solving the linear system 

a 

s ^Rj{X il Hi)H j + R (Xi,fii) = 0, i = l,---,g (2) 



whose solution is 
where 



H = -B~ L V 



H 



H; 



B 



/ Ri(Xi,m) 



Ri(Xi, ^ 



^s(Ai,Mi)\ 



V 



\HgJ \Rl(Xg,Hg) ■■■ Rg(Xg,Hg)J 

Here, of course, we assume that generically det B^O. 



/i*o(Ai,/xi)\ 



Ro(h, Hi) 



\R (Xg,Hg)/ 



(3) 



Theorem 1 Suppose that the variables (Xi,Hi) are separated i.e. they Poisson commute for i ^ j: 

{Xi, Xj} = 0, {Hi, Hj} = 0, {Xi, Hj} = p(Xi,Hi)d~ij (4) 
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Then the Hamiltonians Hi, i = 1 • • • g, defined by eq. (Qj Poisson commute 

{H i ,H j } = 



Proof . Let us compute 

B 1 B 2 {(B- 1 V) 1 ,(B- 1 V) 2 } = {B 1 ,B 2 }(B- 1 V) 1 (B- 1 V) 2 

-{B 1 ,V 2 }(B- 1 V) 1 -{V 1 ,B 2 }(B- 1 V) 2 + {V 1 ,V 2 } 

Taking the matrix element i,j of this expression, we get 

( J B 1 B 2 {(B-V) 1 ,(i?- 1 y) 2 }) = di^iB^Bu^B-'VUB-'V^ 

k,l 

-Sij Y,i B ik, Vi}{B- l V) k - Sij ^{V^BujiB^V)! + SijiVi, F;} = 

k I 

where <5jj occurs because the variables are separated. 



It can hardly be simpler. The only thing we use is that the Poisson bracket vanishes between 
different lines of the matrices, and then the antisymmetry. We did not even need to specify the 
Poisson bracket between Aj and /Uj . The Hamiltonian are in involution whatever this Poisson bracket 
is. This is the root of the multihamiltonian structure of integrable systems. 

Can we make it quantum? Let us consider a set of separated variables 

[Aj, Xj] = 0, fij] = 0, [Xi,/J.j] = p(\i, m)5ij 
We want Baxter's equation, so we start from the linear system 

Y J Rj{\,Vi)H j + Ro{\ i ,iJLi) = Q (5) 

3 

Here the Hj are on the right, and in Rj(\i, fii), Ro(\i,/j>i), we assume some order between A,,^j, 
but the coefficients in these functions are non dynamical. Hence we start from the linear system 

BH = -V (6) 

We notice that we can define unambiguously the left inverse of B. First, the determinant D of B is 
well defined because it never involves a product of elements on the same line. The same is true for 
the cofactor Ajj of the element By (we include the sign (— l) 1 " 1 "- 7 in the definition of A^). Define 

= {B-% = D- l A Jt 

We have 

(B~ l B)ij = 'y^^D~ 1 /± ki B k j 

k 

But Afcj does not contain any element B k i, hence the product A k iB k j is commutative, and the 
usual construction of the inverse of B is still valid. Since the left and right inverse coincide in an 
associative algebra with unit, we have the identities 

{BB~ l )ij = y~] B ik B^l = '^2B ik D~ 1 Aj k = 5ij (7) 

k k 
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We write the solution of eq.@ as 

H = -B~ l V (8) 

Theorem 2 The quantities Hi defined by eq. (^), which solve Baxter's equations eqs.^), are all 
commuting 

II- II, 

Proof . Using that Vk and Vi commute, Vj\ = 0, we compute 

[Hi, Hj] = ^[B-^B/V] (9) 

k,l 



Using 
so that 



k.i 

[A~ l ,B~ x ] = A- 1 B- l [A,B]B- l A- 1 = B^A^^BjA^B- 1 

[ B ik' B ]l l \ = E B ir lB ^[ B rs,B r > s >]B^lB- k } 
rs,r's' 

= E B j>' jB '> 1 [ jBrs ' jBr ' s '] jB ^ ljB s~i 1 



the first term can be written 



it,/ 

= E 2 B ^' B ^ 1 ^ BrS: Br ' S '^ (. B sk B 7n + B sl lB 7'k) V k V l 

Using that [S rs , £? r v] = 5 rr / [B rs , B rs i] and is therefore antisymmetric in ss', and setting 

K ss > = E ( B s>l B 7k + B 7'k B sl l ~ B sl lB s'k - B sk B s'f) V k V l 
k,l 

we get 



fc,Z rss' 

rss' 

= Eg ^ 1 ' B Jr ] [-^rs ! B rs') K ss > 
rss' 

The last two terms in eq.(||) are simpler, we get 

fc,i rsfc 
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The quantities Hi will commute if 



[B-\B~ 1 }=0, Vi,j,r 



(10) 



This is true as shown in the next Lemma. 



The condition eq.(10) says that the elements on the same column of B 1 commute among them- 
selves. In a sense this is a condition dual to the one on B. It is true semiclassically because 

W 1 ' B jr) = Yl B ia B ja'{ B ab, B a>v} B fr B vl = Yl B ia B jai B ab, B ab ,} B^ 1 B^ = 
a,a',b,b' a,b,b' 

where in the last step we use the antisymmetry of the Poisson bracket. We show that it is also true 
quantum mechanically^]. 



Lemma 1 Let B be a matrix whose elements commute if they do not belong to the same line 

Then the inverse B^ 1 of B is defined without ambiguity and moreover elements on a same column 
of B~ l commute 

[Br\Br 1 ] = 



Proof. We want to show that 



denote by pP the vector with components B^, k ^ r. Then we have (with j > i) 

-l) r+i pP A 4 r) A • • • ffi A - - - fiP A • • • pW B ji 

_±y+i+9-j+l R (r) A R {r) 



ri jr r 3 it 



Ki B Jr 



kr 



ft* A (%> A • • • ft) A • • • /J« A • • • PP A ?P B 

-iy+i+9-J+ipp A pP a • • • pP a • • • PP a • • • pp a pPb~} 

-iy+ipP a a • • • pp a • • • pp a • • • pp B ^ 



A ri B. 



r 3 ir 



where in the third line we used eq. ([?]). In the above manipulations, we never have two operators Bij 
on the same line so we can use the usual properties of the wedge product. Moreover it is important 
that the line r is absent in the definition of P^> . Remark that this equation can also be written 
A r iD~ 1 A r j = A r jD~ 1 A r i which is a Yang-Baxter type equation. ■ 



With this Lemma, we have completed the proof of our theorem. It is remarkable that, again, only 
the separated nature of the variables Aj,/Xj is used in this construction, but the precise commu- 
tation relations between Aj,/ij does not even need to be specified. This is the origin of the multi 
Hamiltonian structure of integrable systems, here extended to the quantum domain. 

3 The same Lemma already appeared in jn]]. Our two proofs are different and independent however. We thank B. 
Enriquez for drawing our attention to his work. 
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3 Choosing the right number of dynamical moduli. 



Let us explain how one can set up things in order that the number of dynamical moduli is equal to 
the genus of the Riemann surface. To understand the origin of the conditions we will write, let us 
explain first what happens in the setting of general rational Lax matrices described in |4], 0]. Quite 
generally, a Lax matrix L(X) depending rationally on a spectral parameter A, with poles at points 
Xk can be written as 

L(A) = L + ^L fc (A) (11) 

k 

where Lq = Diag(ai, • • • , ajv) is a constant diagonal matrix and Lk{X) is the polar part of L(X) 
at Xk, ie. Lk(X) = J2r=-n k Lk,r(X — Xk) T . In order to have a good phase space to work with, we 
assume that Lk(X) lives in a coadjoint orbit of the group of N x N matrix regular in the vicinity 
of A = Xk, i.e. 

Lk = (gk^-kgk 1 



Here Ak(X) is a diagonal matrix with a pole of order rik at A = Xk, and gk has a regular expansion 
at A = Xk- The notation ()_ means taking the singular part at A = Xk- This singular part only 
depends on the singular part (Ak)~ and the first nk coefficients of the expansion of gk in powers 
of (A — Xk)- The matrix (Ak)- is an orbit invariant which specifies the coadjoint orbit, and is 
not a dynamical variable. It is in the center of the Kirillov bracket which as shown in Q induces 
the Poisson bracket eq.([|), with p(Xi, Hi) = 1, on the separated variables. The physical degrees of 
freedom are contained in the first n/% coefficients of gk{X). Note however that since Ak commutes 
with diagonal matrices one has to take the quotient by gk — > gkdk where dk(X) is a regular diagonal 
matrix, in order to correctly describe the dynamical variables on the orbit. The dimension of the 
orbit of Lk is thus N(N — l)nk so that L(X) depends on ^ fc N(N — l)rik degrees of freedom. 
Finally, the form and analyticity properties of L(X) are invariant under conjugation by constant 
matrices. To preserve the normalization, Lq, at oo these matrices have to be diagonal (if all the 
Oj's are different). Generically, these transformations reduce the dimension of the phase space by 
2(iV - 1), yielding: 

dimTW = (N 2 - N)^2n k - 2(N - 1) 

k 

The spectral curve is 

JV-l 

r : J R(A,/i) = det(L(A)-/xl) = (- / x) Af +^r (? (AK = (12) 

9=0 

The coefficients r q (X) are polynomials in the matrix elements of L(X) and therefore have poles at 
Xk- The curve is naturally presented as a A^-sheeted covering of the A-plane. We call fJ-j(X) the N 
branches over A. Using the Riemann-Hurwitz formula, we can compute the genus of T 

9 = -^-^ J ~2^n k -N + 1 = -dimM 



It is important to observe that the genus is half the dimension of phase space. So the number of 



action variables occurring as independent parameters in the eq.([12|) should also be equal to g. Let 
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us verify it. Since rj(X) is the symmetrical function <Tj((ii, ■ ■ ■ ,^n), it is a rational function of A. 
It has a pole of order jrik at A = Afc. Its value at A = oo is known since fJ,j(X) — > aj. Hence it can 
be expressed on j parameters namely the coefficients of all these poles. Altogether we have 

\N(N + 1) Y2k n k parameters. They are not all independent however. Above A = Afc the various 
branches can be written: 

^ ( A ) = Yl (X -\ k )n + r6gUlar ( 13 ) 

71=1 

(i) (i) 

where all the coefficients c\ , • • • , Cn fe are fixed and non-dynamical because they are the matrix 
elements of the diagonal matrices (Aj~)-, while the regular part is dynamical. This implies on rj(X) 
that the coefficients of its n& highest order pole terms are fixed. Summing over j, we get Nn^ 
constraints and we are left with ^N(N — 1) parameters, that is g + N — 1 parameters. It 

remains to take the quotient by the action of constant diagonal matrices. The generators of this 
action are the Hamiltonians H n = (l/n)ies\ =ao Ti (L n (\)) dX, i.e. the term in 1/A in Tr(L n (A)). 
Setting 

H {\) = a 3 + b j- + --- (14) 

around the point Qj = (oo,aj), we have H n = X^j ^" 1 ^'- After Hamiltonian reduction these 
quantities are to be set to fixed (non-dynamical) values. So, both aj (by definition) and hi are 
non dynamical. On the functions rj(X) this implies that their expansion at infinity starts as 

(0) (o) (—1) 

rj(X) = Tj + -^j !"•••) with r - and non dynamical. Hence when the system is properly 

reduced we are left with exactly g action variables. 



The constraints eqs.(|l3|, 14) can be summarized in a very elegant way [|. Introduce the differen- 
tial 5 with respect to the dynamical moduli. Then our constraints mean that the differential 5fidX 
is regular everywhere on the spectral curve because the coefficients of the various poles being non 
dynamical, they are killed by 5: 

5(i dX = holomorphic 

Since the space of holomorphic differentials is of dimension g, the right hand side of the above 
equation is spanned by g parameters which are the g independent action variables we were looking 
for. Notice that these action variables are coefficients in the pole expansions of the functions Tj(X), 
and thus appear linearly in the equation of T. Hence eq.(|l2|) can be written in the form eq. ([!]). 
Clearly, these considerations can be adapted by considering more general conditions such as 

— = holomorphic 

fi n X m 



4 Examples. 



Let us show how well known models fit into our scheme. For the hyperelliptic ones, things are so 
simple that we can directly check the commutation of the Hamiltonians. 
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4.1 Neumann model. 



The spectral curve can be written in the form 



2 _ n^ 1 (A-^) _nA) 



^ S 1 : as) 



Performing the birational transformation s = fiQ(X), we get: 

s 2 = Q(A)P(A) (16) 

which is an hyperelliptic curve of genus g = N — 1. The polynomial Q(X) is non dynamical. We 
have (N — 1) independent dynamical quantities, namely the (iV — 1) symmetrical functions of the 
bi, coefficients of P(X). We have 

5u dX = dX = — ^—dX = holomorphic 

P 2fiQ(X) 2s 1 

Asking that a curve of the form eq.(|l^) passes through the g points (Aj,/Xj) determines the poly- 
nomial -P(A). The solution of Baxter's equations 

P(Xi) = Q{h)i4 

simplifies in this case because the matrix B depends only on the Aj. It is equivalent to Lagrange 
interpolation formula: 

P(A) = P(°)(A) + P( 2 )(A) 

with 

P(°)(A) = n(A - A,), P^(A) =£SiW0(A,g m = " 

j llfe^(Ai-A fc J 

Introducing the canonical commutation relations 

[Hj,X k ] = ihSjk 

so that 

[ H , /(A,)] = ihd x J(Xj), [ M 2 , /(X,)} = 2ihd Xj f(X J ) H + (zhfdljiXj) 

We can check that [P(A),P(A')] = is a consequence of 9^.p(°)(A) = — Tik^ji^j ~ ^k)Sj(X), and 
the identities 

SjW^SiW - SjiX'^S^X) = 0, Vn > 

These identities follow from the remark that, if we define the translation operators tjXi = Xj + cr5ij, 
then 

SjWtjSiM - Sj(X')tjSi(X) = ^'^~ X ^ l ^~^ ^ - Ai)(A - A') (17) 

llfc?y(Aj - XkjU^Xi - X k ) 

is independent of a. 
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4.2 Toda Chain. 



The spectral curve can be written in the form |J: 

/i + /i" 1 = 2P(A) (18) 

where 2P(A) = A" +1 — ^7=1 P«A n + • • • is a polynomial of degree (n + 1). The spectral curve is 
hyperelliptic since it can be written as 

s 2 = p 2 (\) - l, with s = n - P(A) (19) 

The polynomial P 2 (A) is of degree 2(n + 1) so the genus of the curve is g = n. The number of 
dynamical moduli is g = n in the center of mass frame YH=i Pi = 0- We have 

— dX = ^--^rdX = — ^-dX = holomorphic 

fl fl — fl 1 s 



Asking that the curve eq.(18) passes through the n points (Aj,//j), we get Baxter's equations. 

2P(A,) = m + v; 1 
Their solution is again given by Lagrange interpolation formula: 

2P(A) = P(°)(A) + PW(A) 

where 

n 

P(°)(A) = (A + £ A,) [](A - A;), P (1) (A) = £ Si(X)(fii + ft 1 ) 

i i=l i 

The polynomial P(°'(A) is of degree n + 1, vanishes for A = Aj and has no A n term. Let the 
commutation relations of the separated variables be given by: 

= q\j/J,j, /UjAj = XiHj, i ^ j 

Then again [P(A),P(A')] = as a result of eq. (|l7|) , where tj is interpreted as tjXj = qXj, and the 
facts that 



S , i (A)i J ± 1 p(°)(A / ) - Sj(X')tf l pW (A) = P(°)(A')5 3 -(A) - p( )(A)5,-(A') 

nw*-*>rw*'-*»> (A + _ 



4.3 A non— hyperelliptic model. 

We consider the model studied in |2|, ||, f?j . The spectral curve can be written in the form: 

R(X,fi) = fi N + t^(X) f i N - 1 + ---t^ N \X) =0 (20) 
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The polynomials fW(A) are such that degree #)(A) < kn - 1 and degree t^ N \X) = Nn - 1 for 
some integer n. The genus of this curve is 

g = l(N-l)(Nn-2) 

Assuming that there is no singular point at finite distance, the homomorphic differentials are 

l \ k 

u) ki = _ ~; , dX, < Z < N - 1, 0<k<(N-l - l)n - 1 

We have 

dA _ ^ 6tW(\)fjL N - k d\ 

This will be holomorphic if (^^(A) = and 

5i«(A) = 5H^\ + ■■■ + SH^i^X^- 1 ^- 1 

Baxter's equations and the commutation of the Hamiltonians where proved in this case, starting 
from the definition of the quantum model through it Lax matrix and the Yang-Baxter equation. 
Our approach gives a very simple proof of this result. 



5 Conclusion. 



We have shown that starting from the separated variables, one can give an easy definition of 
a quantum integrable system. The next step is to reconstruct the Lax matrix and the original 
dynamical variables of the model. While this is a well understood problem in the classical theory 
]l], B], (it is the essence of the classical inverse scattering method), its quantum counterpart will 
require a deeper understanding of the quantum affine Jacobian flC)[| , 

Acknowledgements. We thank D. Bernard and F. Smirnov for discussions. 
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